hep-ph/0612119 
KA-TP-12-2006 



Next-to-leading order QCD corrections to slepton pair 
production via vector-boson fusion 



Partha Konar and Dieter Zeppenfeld 

Institut fur Theoretische Physik, Universitat Karlsruhe, 
P.O.Box 6980, 76128 Karlsruhe, Germany 



Abstract 

Slepton pairs can be produced in vector-boson fusion processes at hadron colliders. 
The next-to-leading order QCD corrections to the electroweak production cross section 
for pp — > £ + £~ + 2 jets at order a s a 4 have been calculated and implemented in a NLO 
parton-level Monte Carlo program. Numerical results are presented for the CERN 
Large Hadron Collider 
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Introduction : 

Among the primary goals of the CERN Large Hadron Collider (LHC) are the discovery of 
the Higgs boson, thus shedding light on the yet unexplored mechanism of electroweak (EW) 
symmetry breaking, and the search for physics beyond the standard model. Within the area 
of Higgs boson studies, vector-boson fusion (VBF) processes have emerged as being highly 
promising for revealing information on the symmetry breaking sector [1]. The prototypical 
process is qq — > qqH, which proceeds via t-channel W or Z exchange. The two scattered 
quarks emerge as forward and backward jets (called tagging jets) which provide a character- 
istic signature for VBF and allow to significantly suppress backgrounds. As a result, VBF 
searches are expected to lead to quite clean Higgs boson signals. 

A natural question is whether vector boson fusion is a useful tool also for the study 
of other signals of new physics. Some recent work has indicated the effectiveness of VBF 
channels in the context of new physics searches, particularly for new particles that do not 
interact strongly. Perhaps the best example [2] is afforded by supersymmetric theories, 
wherein conventional search strategies for neutralinos and charginos may run into difficulty 
at the LHC, for a significant part of the parameter space. The possibility of a slepton search 
has been studied for vector-boson fusion as well [3]. A more recent study [4] on VBF slep- 
ton production using Smadgraph arrived at a substantially smaller cross-section, however, 
which is partly caused by large cancellations among VBF-type diagrams and bremsstrahlung 
diagrams at the Born level. 

The discrepancies between these previous results lead us to a recalculation of the slepton 
pair-production cross section in VBF. The relevant Feynman graphs for this process are 
depicted in Fig. 1 for the tree level contributions. In this approximation, we confirm the new 
results of Ref. [4]. In addition, we also perform a calculation of the NLO QCD corrections 
to this VBF process. The NLO calculation closely follows previous calculations for Hjj and 
Zjj production in VBF in Refs. [5, 6]. It uses the Catani-Seymour subtraction scheme [7] for 
implementing the real and virtual NLO contributions in the form of a fully flexible parton 
level Monte Carlo program. 
Calculation : 

The Feynman graphs contributing to pp — > + 2 jets at tree level are indicated in 
Fig. 1. Considering the possible choices of external quarks or anti-quarks, the sub-processes 
can be grouped into neutral-current (NC) processes, like uc — » uc£ + £~ , and charged-current 
(CC) processes, like us — > dci + i~. For the purpose of calculating the virtual QCD correc- 
tions, the Feynman graphs are divided into Compton scattering type graphs, as in Fig. 1(a), 
and the VBF type graphs as in Fig. l(b-e). The first class (Fig. 1(a) and three additional 
bremsstrahlung diagrams, with the vector boson radiated at the position of the blobs) cor- 
responds to the emission of the external vector boson from one of the two quark lines. The 
VBF graphs represent VV — > £ + £~. Here, V stands for a t-channel 7, Z or W boson. For 
selectron or smuon production one expects a negligible contribution from Fig. 1(b). We 
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(a) 



(b) 




Figure 1: Generic LO parton level diagrams leading to slepton pair-production through elec- 
troweak VBF at hadron colliders. 

do include this Higgs exchange contribution for stau pair production, however, anticipating 
strong enhancements of the couplings to the Higgs bosons at large tan /3. 

Contributions from anti-quark initiated t-channel processes such as uc — > uc£ + £~, which 
emerge from crossing the above processes, are fully taken into account. On the other hand, 
two additional classes of diagrams which can appear in case of identical quark flavors, are 
simplified in our calculation. The first concern s-channel exchange diagrams, where both 
virtual vector bosons are time-like. These diagrams correspond to vector boson pair produc- 
tion with subsequent decay of one neutral vector boson to £ + £~ while the other one decays 
into a quark-antiquark pair. These contributions can be safely neglected in the phase-space 
region where VBF can be observed experimentally, with widely-separated quark jets of large 
invariant mass. The second class corresponds to w-channel exchange diagrams which are 
obtained by the interchange of identical final state (anti) quarks. Their interference with the 
t-channel diagrams is strongly suppressed for typical VBF cuts and therefore neglected in 
our calculation. Color suppression further reduces any interference terms. 

Throughout our calculation, fermion masses are set to zero and external b- and t-quark 
contributions are neglected. For the Cabibbo-Kobayashi-Maskawa matrix Vckmi we use a 
diagonal form equal to the identity matrix. This yields the same results as a calculation 
using the exact Vqkm when the summation over all quark flavors is performed. 
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The computation of NLO corrections is performed in complete analogy to Ref. [6]. For the 
real-emission contributions, we consider the diagrams with a final-state gluon by attaching 
the gluon to the quark lines in all possible ways. As a result one obtains two distinct, non- 
interfering color structures which correspond to gluon emission off the upper or off the lower 
quark line in Fig. 1. Subprocesses with an initial gluon are obtained by crossing the final 
state gluon on a given quark line with the incident quark or anti-quark of this same quark 
line. As a result only one color structure exists for initial gluons. The other color structure 
would correspond to an s-channel process of the type gq — > VVq, which has been neglected 
also at Born level. 

All amplitudes are evaluated numerically using the amplitude techniques of Ref. [8]. The 
calculation is simplified by introducing the leptonic tensors T v and Ly V , which describe the 
effective polarization vector of the final state decay V(q) — > £~ (p-\)i + (p-i) , 

V£ 

r£(pi,p 2 ) = , - V2 9r 2 — = (Pi-P 2 ) Q , (i) 

(pi + P2) 2 -m v + imv^v 

and the effective sub-amplitude for the process VfV^ — > The leptonic tensor T v is 

common to real emission graphs and Born graphs appearing in the Catani-Seymour subtrac- 
tion terms and needs to be calculated only once at a given phase space point, independent of 
the crossing of the colored partons. Similarly, Lyy is only needed for two distinct momentum 
flows (gluon attached to the upper or to the lower quark line) at any phase space point. It is 
calculated in the complex-mass scheme [9] which implements the Breit-Wigner propagators 
of the resonant Z-boson in a gauge invariant way. 

At NLO, we have to deal with singularities in the soft and collinear regions of phase space 
which are regularized in the dimensional-reduction scheme [10] with space-time dimension 
d = 4 — 2e. The cancellation of these divergences with the respective poles from the virtual 
contributions is performed by introducing the counter terms of the dipole subtraction method 
[7]. Since these divergences only depend on the color structure of the external partons, the 
analytical form of subtraction terms and finite collinear pieces encountered for VBF slepton 
pair production, in terms of the respective Born amplitude, is identical to the ones given in 
Ref. [5]. 

The virtual corrections to the amplitudes arise from a virtual gluon emitted and re- 
absorbed by either the upper fermion line or by the lower fermion line. For both contributions 
the resulting virtual amplitude, M.y, can be expressed in term of a divergent part, which is 
fully factorisable in terms of the original Born amplitude, A4b, and a finite part, A4y, 

My = M B ^>C F (^)'r(l +«)(—-! + C*<) + My . (2) 

Here, the first term gets contributions from virtual QCD corrections to all types of Feynman 
graphs as in Fig. 1 but the finite second part originates from virtual QCD corrections to 
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Figure 2: Scale dependence of the pp — > £ + £ jjV cross section at the LHC for (a) left-type 
slepton and (b) right-type slepton at NLO and LO with the cuts of Eqns.(3, 4)- The slepton 

= 200 GeV. 



masses are m^ L = m 



only those Feynman graphs where two electroweak bosons are attached to the same fermion 
line, as shown in Fig. 1(a). The full expression of this finite part can be expressed in terms 
of the finite parts of the Passarino-Veltman [11] Bq, Co and Dij functions and was given in 
Eq. (Al) of Ref. [6] for the analogous case of V — > decay: simply replace one of the two 
polarization vectors e° of Eq. (Al) by the slepton current Ty of Eq. (1). 

The results obtained for the Born amplitude, the real emission and the virtual correc- 
tions have been tested extensively. For the tree- level amplitudes (Born and real emission), 
we have performed a comparison to the fully automatically generated results provided by 
Smadgraph [4] and confirmed their equality numerically. We also checked the invariance of 
the Born cross-section under Lorentz transformations. Furthermore, gauge invariance has 
been confirmed for the external gluon, within the numerical accuracy of the program. 
Results and Discussions : 

The cross-section contributions discussed in the previous section are implemented in a 
fully-flexible parton-level Monte Carlo program for EW £ + £~jj production at NLO QCD 
accuracy. The program is very similar to the ones for Hjj, Vjj and VVjj production in 
VBF described in Refs. [5], [6] and [12]. We use the CTEQ6M parton distributions with 
a s (mz) = 0.118 at NLO, and CTEQ6L1 distributions for all LO cross sections. We chose 
m z = 91.188 GeV, m w = 80.423 GeV and G F = 1.166 x 10~ 5 GeV" 2 as electroweak 
input parameters. Thereof, oqed = 

1/132.54 and sin 2 W = 0.22217 are computed via LO 
electroweak relations. To reconstruct jets from final-state partons, the kT algorithm is used 
with resolution parameter D = 0.8 [13]. Throughout, we assume a pure decay of the 
sleptons, whenever decay distributions are being discussed. 
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Figure 3: Distributions for (a) daughter lepton invariant mass M\\ and (b) missing transverse 
momentum pT m i SS at NLO (solid red) and LO (dashed black). Left type slepton production 
in the SPS la scenario (mi = 202 GeV) is considered. Renormalization and factorization 
scales are taken as hr = fip = 1v ■ 



Partonic cross sections are calculated for events with at least two hard jets, which are 
required to have 

p Tj >20GeV, |%|<4.5 (3) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four-vector sum of massless partons of pseudo-rapidity \rj\ < 5. These cuts ensure a finite LO 
differential cross section for £ + £~jj production, since they enforce finite scattering angles for 
the two quark jets. The two reconstructed jets of highest transverse momentum are called 
'tagging jets'. At LO, they are the final-state quarks which are characteristic of vector-boson 
fusion processes. Backgrounds to VBF are significantly reduced by requiring a large rapidity 
separation of the two tagging jets. We therefore impose the cut 

&y jj = \y jl -y j2 \>±2 (4) 

Within the above cuts we have calculated the £ + £~jj cross sections at LO and at NLO 
for the SPS la parameter point where slepton masses are given by mi L = 202 GeV, ^i R = 
144 GeV. This point can be parameterized by the mSUGRA model with m = 100 GeV, 
mu2 = 250 GeV, A = —100 GeV, tan/5 = 10 and positive \i [14]. We find production cross 
sections of 0.0536 (0.0532) fb for £ L production and 0.0242 (0.0249) fb for £ R production 
at NLO (LO) when setting renormalization and factorization scales to hr = \ip = qv- 
Unfortunately, expected cross sections at the LHC are quite small in general, not exceeding 
0.1 fb for slepton masses above 150 GeV for left-handed sleptons and for slepton masses 
above 80 GeV for right-handed sleptons. In order to compare £l and Ir cross sections more 
directly, we have calculated their total production cross sections, within the above cuts, for a 
mass of 200 GeV in both cases. Fig. 2 illustrates the dependence of these total cross sections 
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Figure 4: Distributions of the azimuthal angle between the two tagging jets , <f>jj, at NLO 
(solid red) and LO (dashed black). Left type slepton production in SPS la scenario (W| = 
202 GeV) is considered. Renormalization and factorization scales are taken as jiR = hf = 

qv- 

on the renormalization and factorization scales, hr and /xp, which are taken as multiples 
of the momentum transfers, qv, of the t-channel electroweak bosons in Fig. 1, [j,r = £r<jv, 
fip = £fQv- This choice takes into account that at both LO and NLO the VBF process 
can be viewed as a double deep inelastic scattering event, for which the momentum transfer 
carried by the exchanged electroweak boson is a natural scale choice. It leads to K-factors 
close to unity for both total cross sections and distributions. 

The LO cross section, o^o, only depend on \ip = ^Fqv- By varying the scale factor 
^f = £ in the range 0.1 - 10, the value of o~lo changes by around a factor of two, indicating 
a substantial theoretical uncertainty of the LO calculation. The strong scale dependence is 
reduced at NLO. For a^Lo, we show three different cases: £p = (,r = £ (solid red line), 
£f = £, £r = 1 (dot-dashed blue line), and £f = 1, £r = £ (dashed green line). The latter 
curve illustrates clearly the weak dependence of ct^lo on the renormalization scale, which 
can be understood from the fact that a s (fiR) enters only at NLO. Also the factorization- 
scale dependence of the full cross section is low. In our following study we fix the scales at 
A*f = — Qv, unless noted otherwise. 

Two examples for distributions are given in Fig. 3. We show the distributions for (a) the 
invariant mass, M U) of the two charged daughter leptons in the decay £ + £~ — > l + l~PTmi SS i an d 
(b) the missing transverse momentum, PT m iss- Results are shown at both LO and NLO and 
are virtually indistinguishable with the scale choice fip = ^r = qv ■ For the illustration in 
Fig. 3, left type slepton production (pp — ► i^iJ^jjX) in the SPS la scenario (to| = 202 GeV) 
is considered. 

Within the same set of model parameters, the distribution in the azimuthal angle be- 
tween the two tagging jets, fyj, is shown in Fig. 4. One finds a characteristic dip at 90°, 




7 



0.0008 



2 0.0006 



t— 
f 



0.0004 - 




-NLO Ratio :[n F =q v , m z l 
1.2 - LO Ratio :[^p=q v , m z l 



K factor : Hp^Qv 
K factor : |lp=m z 



1.4 



1 



O 

o 



0.8 - 
1-2 0.6 - 



(I 



50 



mo 



150 



200 



50 100 150 

p T -min tags [GeV] 



50 100 

p T -min tags [GeV] 



Figure 5: (a) Distribution of the minimum transverse momentum of the two tagging jets, 
PT,min(tags) , at NLO (solid red) and LO (dashed black) for the scale choice Hr = [If = Qy ■ 
(b) K factor as defined in Eq. (5) after fixing the factorization scale at the momentum transfer 
fip = qy (solid red) and at Lip = mz (dashed blue). The inset shows the ratio between these 
two choice of factorization scale at NLO (solid red) and at LO (dashed black). Left type 
slepton production in the SPS la scenario (mj> = 202 GeV) is shown. 



a feature which otherwise is only found in Hjj production when the Higgs boson couples 
to gauge bosons via a HV^V^ operator in the effective Lagrangian [15]. Hjj or Zjj pro- 
duction via VBF in the SM produces a fairly flat <ftjj distribution, while Hjj production 
via gluon fusion exhibits a structure very similar to the one shown in Fig. 4. Since it has 
been suggested that the <pjj distribution in VBF events be used in distinguishing SM Higgs 
couplings from anomalous couplings, the possibility that a dip at 90° might also be produced 
by the production of two charged scalars should be kept in mind, should such a feature be 
discovered at the LHC. 

A distribution which distinguishes slepton pair production from many other VBF pro- 
cesses is the minimum transverse momentum of the two tagging jets, dcr/dpT,min(tags), which 
is depicted in Fig. 5(a). Due to a significant contribution from t-channel photon exchange, 
this distribution falls quite steeply. An even steeper fall-off is found for right type slepton 
production (pp — > i^i^jjX), which can be understood from the fact that £r has no coupling 
to the W ±,Q eigenstates. 

The shape of the above distribution at LO can differ significantly from the respective 
NLO result when scales other than Hr = Hf = Qy are used. This is emphasized in Fig. 5(b), 
where we show the dynamical K factor, defined as 

da^Lo/dx 



K(x 



(5) 



daLo/dx 

for the two choices tip = Qy ari d Hf — m z (and lir = qy for the NLO curves in both cases). 
While the NLO cross sections differ very little when switching between the two scale choices 
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(see inset), the effect on the LO cross sections is quite sizable, approaching a 20% effect at 
p T j ~ 200 GeV. This is to illustrate that the choice fj, F = q v minimizes the NLO corrections 
in most distributions, by producing a LO prediction which is close to the true NLO result. 
Summary and Conclusions: 

In this paper we have presented results for EW slepton pair production at NLO QCD 
accuracy, obtained with a new parton-level Monte Carlo program. The integrated cross 
sections for this process are consistent with the results of Ref. [4] and show a very moder- 
ate K factor. While NLO results are quite stable against scale variations, LO results can 
change substantially. We find that the higher order QCD corrections are minimized by the 
scale choice hf = Qv at LO, where qv is the momentum transfer carried by the t-channel 
electroweak bosons. 

A second observation concerns the distribution of the azimuthal angle separation between 
the two tagging jets in VBF events. The VBF production of two scalars, as considered here, 
produces the same type of dip at 90° as is otherwise observed only for Higgs production with 
loop induced couplings to the fusing vector bosons. 
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